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What can break the Wandzura-Wilczek relation? 
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We analyze the breaking of the Wandzura-Wilczek relation for the <?2 structure function, empha- 
sizing its connection with transverse momentum dependent parton distribution functions. We find 
that the relation is broken by two distinct twist-3 terms, and clarify how these can be separated 
in measurements of double-spin asymmetries in semi-inclusive deep inelastic scattering. Through a 
quantitative analysis of available <?2 data we also show that the breaking of the Wandzura-Wilczek 
relation can be as large as 15-40% of the size of Q2- 



I. INTRODUCTION 

The spin structure of the nucleon remains one of the 
most challenging and controversial problems in hadronic 
physics [2, [2, 3]. Experimentally it is now known, through 
many careful measurements of the nucleon's g\ structure 
function, that quarks carry only some 30% of the pro- 
ton's longitudinal spin, a feature which is now qualita- 
tively understood (4J. Moreover, polarized pp scattering 
observables [B[ and open charm production in deep in- 
elastic scattering [6| suggest that gluons carry an even 
smaller fraction of the longitudinal spin. Presumably, 
the remainder arises from quark and gluon orbital angu- 
lar momentum. 

Although less attention has been paid to it, there are 
a number of intriguing questions associated with the 
transverse spin structure of the nucleon. An example 
is the study of the g 2 structure function, which only 
in recent years has been probed experimentally with 
high precision. Unlike all other inclusive deep-inelastic 
scattering (DIS) observables, the g 2 structure function 
is unique in directly revealing information on the long- 
range quark-gluon correlations in the nucleon. In the 
language of the operator product expansion (OPE) these 
are parametrized through matrix elements of higher twist 
operators, which characterize the strength of nonpertur- 
bative multi-parton interactions. (In the OPE "twist" is 
defined as the mass dimension minus the spin of a local 
operator.) In other inclusive structure functions higher 
twist contributions are suppressed by powers of the four- 
momentum transfer squared Q 2 , whereas in g 2 these ap- 
pear at the same order as the leading twist. 

As discussed by Wandzura and Wilczek [7| , the leading 
twist contribution to the g 2 structure function, which is 
denoted by g^ w , can be expressed in terms of the leading 
twist (LT) part of the g\ structure function, 
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where xb is the Bjorken scaling variable, and we suppress 
the explicit dependence of the structure functions on Q 2 . 
The Wandzura-Wilczek (WW) relation asserts that the 
total g 2 structure function is given by the leading twist 



approximation JT]), 
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which would be valid in the absence of higher twist con- 
tributions. In this case the g 2 structure function would 
also satisfy the Burkhardt-Cottingham (BC) sum rule 
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Its violation would also signal the presence of twist-3 or 
higher contributions. Unlike the WW relation, however, 
the validity of the BC sum rule (which is yet to be con- 
clusively demonstrated experimentally (slUH) would not 
necessarily imply that higher twist terms vanish [ill, [l^] • 
In this paper we explore the physics that can lead to 
the breaking of the WW relation in QCD, preliminary 
results for which have appeared in Ref. [HI]. In Sec. II 
we present a detailed theoretical analysis of quark-quark 
and quark-gluon-quark correlation functions, and discuss 
the so-called Lorentz invariance relations and equations 
of motion relations. From these we show that the WW 
relation is valid if pure twist-3 and quark mass terms 
are neglected, in agreement with OPE results. We find 
that there are two distinct contributions with twist 3, 
denoted by <?t and gx, which correspond to two different 
"projections" of the quark-gluon-quark correlator. An 
explicit demonstration of our findings is made for the 
case of a point-like quark target, which shows that the 
twist-3 terms can in principle be as large as the twist-2 
terms. 

In Sec. Ill we discuss the phenomenology of the WW 
relation for both the proton and neutron, and find that 
the available data from SLAC and Jefferson Lab indicate 
a breaking of the WW relation at the level of 15-40% 
of the size of g 2 within the 1-a confidence level. The 
two twist-3 terms can be separated by measuring, in ad- 
dition to g 2 , the function <?[^ in semi- inclusive DIS, as 
we outline in Sec. IV. There we explain the importance 
of measuring the two twist-3 functions <?t and <?t sepa- 
rately, and the insight which this can bring, for example, 
to understanding the physics of quark-gluon-quark corre- 
lations [HI j or to determining the QCD evolution kernel 
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for <72 and the large momentum tails of transverse mo- 
mentum distributions (TMDs). 

Finally, in Sec. V we briefly summarize our find- 
ings. Some technical details for the analysis with a non- 
lightlikc Wilson line and the model calculation of parton 
correlation functions are presented in the appendices. 



II. THEORETICAL ANALYSIS 

In this section we set forth the framework for our anal- 
ysis of the WW relation by first defining quark-quark 
correlation functions and examining their most general 
Lorentz and Dirac decomposition. This is followed by 
a discussion of quark-gluon-quark correlators, and of 
the Lorentz invariance and equations of motion relations 
from which a generalization of the WW relation is de- 
rived. 



A. Parton correlation functions 

The quark-quark correlator for a quark of momentum 
A: in a nucleon with momentum P and spin S is defined 
as 
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where the quark fields ipf are labeled by the flavor index 
a and Dirac index i. For ease of notation, the Dirac and 
flavor indices will be suppressed in the following. The 
operator W" ^ represents a Wilson line (or gauge link) 
from the origin to infinity along the direction specified 
by the vector v, and is necessary to ensure gauge invari- 
ance of the correlator. The gauge links contain transverse 
pieces at infinity 15l 3 and their precise form depends 
on the process lid]. In a covariant gauge, the de- 
pendence of the correlator <I> on v is evident from the 
presence of the Wilson line in the direction conjugate to 
v. In light-cone gauges the vector v is orthogonal to the 
gauge field A, v-A = 0, and the dependence on v appears 
explicitly only in the gauge field propagators. 

In tree-level analyses of semi-inclusive DIS (SIDIS) 
PJ, [H or the Drell-Yan process [U, [H HI] v is iden- 
tified with the light-cone vector n_, where r?_ = = n\ 
and n_ • ?i+ = 1, with n + the corresponding orthogonal 
light-cone vector proportional to P (up to mass correc- 
tions). However, factorization theorems beyond tree-level 
[13, HE HE H3] demand a slightly non-lightlike vector v 
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in order to regularize light-cone divergences. We leave 
a more detailed discussion of the effect of the choice of 
v to Appendix [A] and consider v = n_ unless otherwise 
specified. 

The correlator $ can be parametrized in terms of struc- 
tures built from the four vectors P, S, k and v. Its full 
decomposition has been studied in Ref. (28| (and further 
generalized in Ref. [29j|). It contains 12 scalar functions 
Aj already known from Refs. [lEIIE], and 20 scalar func- 
tions Bi which are multiplied by factors depending ex- 
plicitly on v, which were first introduced in Ref. 3l| and 
called parton correlation functions (PCFs) in Ref. [27l ]. 
For brevity we consider only those terms of the full de- 
composition [28[ which are necessary for the present anal- 
ysis, 
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where the nucleon mass M is explicitly included to ensure 
that all PCFs have the same mass dimension. (Any other 
hadronic scale, such as Aqcd, can be chosen, but we 
follow the choice used in the TMD literature |19j.) 

The PCFs Ai and £^ are in principle functions of the 
scalar products P • k, k 2 , P-v, k-v and v 2 . However, be- 
cause the correlator $ is invariant under the scale trans- 
formation v — > Aw, where A is a constant, the PCFs can 
only depend on ratios of the scalar products, P • k, k 2 
and k ■ v/P ■ v. We therefore choose the PCFs to depend 
on the parton virtuality r = k 2 , on a = 2P • k, and on 
the parton momentum fraction x — k ■ n~ / P ■ n_ . We 
emphasize that the explicit dependence on x is induced 
in general by the v dependence of the correlator <I>. 

These considerations apply even when the correlator 
is integrated over the parton transverse momentum, and 
in fact the Bi terms give contributions also to standard 
collinear parton distribution functions (PDFs), such as 
the helicity distribution — see Eq. (|2"2")l below. However, 
when the correlator is fully integrated over d A k the Bi no 
longer contribute; indeed 



d A k $(fc, P, S; v) = (P,S\ V(0) Vi(0) | P S) , (6) 



and the dependence of the integral on v disappears be- 
cause W^W^^l. 



In TMD factorization the relevant objects are the integrals of &(k, P, S; v) over k = k^n 1 ^, 
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It is also useful to define the fey-integrated correlators 
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where LC refers to the correlators in the light-cone gauge. The correlator <j>g actually depends on the detailed form of 
the Wilson line, and changes, for example, between the SIDIS and Drell-Yan processes. However, for our discussion 
this will not be relevant and we can consider the average between the correlator for SIDIS and Drell-Yan 151 1. 

I 



For any correlator, we can introduce the Dirac projec- 
tions 



1 



$ [r] = -Tr[r$] , 
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where T is a matrix in Dirac space. The transverse mo- 
mentum dependent parton distribution functions then 
appear as terms of the general decomposition of the pro- 
jections $[ r l(x, fey), the full list of which can be found in 
Refs. [USl]. Usually a TMD is defined to have "twist" 
equal to n if in the expansion of the correlator it appears 
at order (M/P+)™ -2 , where P + = P^nt. In physical 
observables, TMDs of twist n appear with a suppression 
factor (M/Q) n - 2 compared to twist-2 TMDs. We finally 
note that at present TMD factorization for SIDIS has 
been proven for twist-2 TMDs only 24], and problems 
are known to occur at twist 3, indicating that the for- 
malism may not yet be complete (32l l33j . 

For the following discussion we shall need the defini- 
tions of certain TMDs (note that here and in the follow- 
ing a is restricted to be a transverse index) [20| 
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where Sl = S + AI/P + , and the transverse tensors g^f 



and e^ p arc defined as 
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(14) 
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For the fey-integrated distributions, we correspond- 
ingly have 



^ + ^(x) = S L g 1L (x), 
$ ^ Q+ 75]( a ,) =s$h 1 (x), 

*P + T ] (x) = S$Mg$(x), 
where for any TMD / = f(x, fey) we define 
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To avoid confusion with the structure function gi, here 
we use the notation g\L also for the helicity-dependent 
PDF, contrary to what is used in some of the TMD lit- 
erature [20l | . 

The connection between the TMDs and the Ai and Bi 
amplitudes has been worked out in detail in the Appendix 
of Ref. [IIS] for v = n_. In Appendix lAl we extend these 
results to a non-lightlike vector v. We shall not repeat 
here the calculations but only quote the results relevant 
for our discussion, namely 
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x (A 7 + xA s ) , (23) 
gx{x, k T ) = J dadr S(t — xa + x 2 M 2 + k T ) 

x(-Ae ^ As). (24) 

As anticipated, we see that Bi terms appear also in the 
function giL, which survives if the correlator is integrated 
over kx- 

B. Lorentz invariance relations 

From the preceding discussion, using the techniques 
discussed for example in Ref. [3(3] , it is possible to derive 
the so-called Lorentz invariance relation (LIR) 



9t(x) =9xl{x) + -r.9yr( x ) + 9t(x), 
where the function gx is given by 
gx{x) = J d 2 k T dad,TS(T - xcr + x 2 M 2 + k T ) 



(25) 
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The proper operator definition for gx can be traced back 
to Ref. [35[ (see also [H,[13)) and requires the introduc- 
tion of the twist-3 quark-gluon-quark correlator 
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where _F +a is the gluon field strength tensor, k' is the 
gluon momentum, and x' = k' -n^/P-n^. Note that this 
correlator has been discussed in slightly different forms 
in Refs. [H, [H, Hi| H(| , for example. It can be expanded 
in terms of four scalar functions Gp, Gp, Hp and Ep 
according to [H, H(| 



i$ F (x,x') = — 
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Hermiticity and parity invariance impose that these func- 
tions are real and either odd or even under the inter- 
change of x and x 1 |4fJ, 

Gf{x, x') — Gf(x', x) , Gf{x, x) — — Gf(x', x) , (29) 
Ef{x, x') — Ef(x', x) , Hf{x, x') — — Hf(x', x) . (30) 
We can then express the function 7jx as 

,i<S> F [ ^\x',x) 



MS£g T (x) 



dx 



{x-x'Y 



MSZ V / dx' 



G F (x,x')/(x-x') 



(31) 



where V denotes the principal value integral. (The 
need for the principal value was apparently overlooked 
in Refs. (35 , l37l|.) The imaginary part arising from the 
pole at x = x' cannot give a contribution to the LIR in 
Eq. (f25|) , but rather contributes to a LIR involving the 

functions fx and f^f , which wc do not discuss here. 
We note that gx is a "pure twist-3" function, being part 
of the twist-3 correlator of Eq. (|2~7| . Since the integrand 
in Eq. (|3ip is antisymmetric in x <-> x' , one obtains the 
nontrivial property 

f i 

dxg T {x)=Q. (32) 

D 

In some analyses [lo|, EU #t was believed to van- 
ish because (i) the Bi parton correlation functions were 
not taken into account, (ii) the partial derivatives in 
Eq. (j2"B"|) were neglected since an explicit x-dependcncc 
of the PCFs is generated only through the additional v- 
dependence, (hi) the boundary terms like the last terms 
in (j26|) were neglected. However, none of these assump- 
tions is justified, as we show explicitly in a quark-target 
perturbative calculation in Appendix [Bj We can fur- 
ther draw some model-independent conclusions about the 
boundary terms by comparing them with the expres- 
sion for gix in Eq. (I23p . Positivity bounds imply that 
I^tSitI < M\kx\fi |42j], which is sufficient to guarantee 
that the k T = boundary term indeed vanishes. How- 
ever, since gvr behaves as l/k T at large kx [HI, the 
boundary term at k T — oo cannot be neglected. 

If 'gx is nonetheless neglected, it is possible to ex- 
press the twist-3 function gx in terms of the twist-2 
functions <?il and g\x fl9L f30j ] . Relations of this kind 
have been often mistakenly called Lorentz invariance re- 
lations [l9|, l3fj, |43j , but should not be confused with the 
correct Lorentz invariance relations such as in Eq. (|25p . 

In the literature, model calculations have been used 
to argue that the pure twist-3 terms are not necessarily 
small [ill, HI]- For example, gx can be computed per- 
turbatively in the quark-target model of Refs. [13, H]. 
Using Eqs. (38), (40) and (42) of Ref. [13 one finds 

a (I 2 

gx(x) - g 1L (x) = -^C F In ^ [2x - 6(1 - x)] , (33) 
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where Cf = 4/3, /i is an infrared cutoff, and from 
Eq. ([25]) one has 



9t(x) 



2tt 



Cfc In ^[1-5(1 



/' 



(35) 



From this calculation one can see that gj* is comparable 
to the size of the other twist-2 functions. Moreover, its 
lowest moment vanishes, so that the nontrivial require- 
ment of Eq. (|32| is fulfilled. In Appendix [Cl we confirm 
the above result (for x < 1 only) starting directly from 
the definition in Eq. (f3Tj) . 



C. Equations of motion relations 

The equations of motion (EOM) for quarks, Iptp — mip, 
imply further relations between twist-2 and pure twist-3 
functions (namely, between qq and qgq matrix elements) . 
They are referred to as "equations of motion relations" , 
and for the case of interest here are expressed as 

9rr( x ) = xgr(x) - xg T (x) - — h t (x) , (36) 



where 
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The full list of EOM relations can be found in Ref. [20(. 

Using Eq. (|36|) to eliminate g^(x) in Eq. ([25]) . one 
finds the differential equation 

d I „ m h\\ „ m hi ^ 

X Tx \3T-9T-^-'j+9i L -9 T - Ji -+9T = Q. 

(38) 

Assuming that the relevant functions are integrable by 
fx(dy/y) an d solving for gx one finds 

9t{x)= [ —(giL(y) + 9T(y) 
Jx y v 

+ g* T (x) + ^(h 1 /x)*(x) , (39) 



where we have introduced the shorthand notation 
/*(*) = f(x) ~ f *»f(y) = - f 1 *L± [yf{y)\ 

Jx y Jx y dy 



(40) 



Note that if the integrals over x and y can be exchanged, 
the function / satisfies 



dxf*(x) = 



(41) 



In general, however, this is not necessarily true, as 

stressed in Refs. [ill fl2j|. 

In DIS on a quark-target, ot can be computed using 
Eqs. (38) and (43) of Ref. $7$, giving 
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Again we see that the twist-3 function gx has a size com- 
parable to that of the other twist-2 functions. 



D. Breaking of the Wandzura— Wilczek relation 

The hadronic tensor relevant for spin-dependent DIS 
structure functions is given by the standard Lorentz de- 
composition 

1 



W^(P,q) = —e^q p 



P-q 

Sagi(xB-,Q 2 ) 



P-q / 



(44) 



where q is the momentum of the exchanged photon and 
xb = Q 2 /(2P ■ q) is the Bjorken variable. In general the 
structure functions g\ and gi in Eq. (|44[) are functions 
of the physical (external) variable xb and are given by 
convolutions of the hard 7*-parton scattering coefficient 
functions and the relevant PDFs. At leading order in a s , 
and including terms up to twist 3, they can be expressed 
in terms of the distributions g\ L and g^ (where we now 
explicitly include the flavor index a) introduced above as 

9i(x) = lj2 e «9 a iL(x) , (45) 

a 

9l (x)+g2(x) = -Y,elg£(x), (46) 

a 

where for simplicity we have suppressed the Q 2 depen- 
dence. This then enables the difference between the full 
172 structure function and the WW approximation |T]) to 
be written as 

g 2 (x) - gf w (x) 

= \ £4(Vr(z) + Yi [K,xy{x)+ f x j rTiy) ) ' 

(47) 

which represents the breaking of the WW relation. Note 
that the right-hand-side of Eq. (f4"T|) contains a quark mass 
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term and two pure twist-3 terms. This is the main result 
of our analysis. 

From Eq. (|4T|) the x integral of the pure twist-3 func- 
tion containing and the mass term vanish. Using 
Eq. (|32"j) , and assuming that g%, is regular enough to ex- 
change the x and y integrals, we see that the ~g^ term 
also vanishes. This implies that the above expression for 
(72 satisfies the Burkhardt-Cottingham sum rule, Eq. 
which is not in general guaranteed in the OPE [IJ, [lj 

To obtain the WW relation one must neglect quark 
mass terms compared to the hadron mass (which can 
be reasonably done for light quarks), and either neglect 
both of the pure twist-3 terms, or assume that they can- 
cel each other. The explicit quark-target perturbative 
calculations show that such a cancellation does not take 
place in general, and that the size of the WW breaking 
term can be comparable to the size of , 



g? w (x) = l-5(l-x) 
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x 2 



2x 2 



(48) 



g 2 (x)-g^ w (x) =5(l-x)-l+ ^C F In % 

2n fjb- 



log 



[l-xf 



+ -S(l-x) + 



(49) 



To obtain the above expressions we again made use of the 
results in Ref. [37j]. Note that both gg™ and the total 
<?2 structure function in the quark-target model respect 
the BC sum rule. 



III. CONSTRAINTS FROM DATA 

It is often stated in the literature (see e.g. Ref. (45| ) 
that the WW relation holds experimentally to a good 
accuracy. While there are certainly indications that this 
may indeed be so 0, [Toj , it is important to quantify the 
degree to which this relation holds and place limits on 
the size of its violation. This is the focus of this section. 

We define the experimental WW breaking term 
A ex (xB) as the difference between the experimental data 
and ff2 ww , 



A C x(x B ,Q 2 )=g c 2 *(x B ,Q 2 



WW/ 



2 ) 



(50) 



with the Wandzura-Wilczek term computed using the 
LSS2006 (set 1) fit of the g x structure function [46]. The 
fit was performed including a phenomenological higher- 
twist term and target mass corrections in order to extract 
the pure twist- 2 contribution, g\ T '. Using parametriza- 
tions of Q\ w hich do not account for the 1 /Q power cor- 
rections 47, 48] would risk inadvertantly including spuri- 
ous higher twist contributions when computing the WW 



approximation. We will demonstrate the impact of this 
difference by comparing our g^ with (g 2 A ' W )' computed 
using the total g\ instead of g\ T in Eq. fTj. 

For proton targets we consider data from the SLAC 
E142 |4j| and E155x Q experiments, while for the neu- 
tron only the high-precision data sets from the SLAC 
E155x [|, and Jefferson Lab E99-117 and EOl-012 
[5~fl ] experiments, obtained using 2 H or 3 He targets, are 
included. We checked explicitly that including the lower- 
precision data sets from Refs. (491 . l52l . [53j does not alter 
the fit results, except for artificially lowering the x 2 val- 
ues due to the much larger errors compared to the higher- 
precision data sets. In total, there are 52 data points for 
the proton and 18 points for the neutron, which are used 
separately to fit the WW breaking term A. Systematic 
errors, when quoted, were added in quadrature. For the 
shape of A we choose the form 

A(x B ,a, 13) = a(l - x B f ((/3 + 2)x B - l) , (51) 

which vanishes at ig = 1, has no divergences at xb = 0, 
fulfills the BC sum rule, and only has a single node. We 
do not consider its Q 2 QCD evolution. The evolution of 
gi has been studied numerically in Ref. [54| in the limit 
of a large number of colors. Most of the data considered 
lie in the range 1 GeV 2 < Q 2 < 10 GeV 2 where the effect 
of QCD evolution is rather mild, as indicated also by the 
results of the EOl-012 experiment [5l| . 

The goodness of the fit is estimated using the \ 2 func- 
tion 



„ \A(x B i) 



{XBi)\ 



(52) 



To quantify the size of the breaking term A compared to 
g^ w we define, for any interval [x^ m , x^ ax ], the ratio of 
their quadratic integrals 



dyx 2 B A 2 {x B ) 
dyx%gl(x B ) 



(53) 



with y = log(xs). The value of r is a good indicator of 
the relative magnitude of A and g2, which change sign 
as a function of Xb- In practice we compute r at the 
average kinematics of the E155 experiment @. For the 
proton, we consider three intervals: the entire measured 
xb range, [0.02,1]; the low-^s region, [0.02,0.15]; and 
the high-xs region, [0.15,1]. For the neutron, due to 
the limited statistical significance of the Iow-ccb data, we 
limit ourselves to quoting the value of r for the large- a; b 
region, [0.15,1]. 

The results of the fits are presented in Table U and Fig- 
ure [T] The proton fit displays a positive WW breaking at 
large-xs and a negative breaking at small-x^. The size 
of the breaking term is typically 15-35% of the size of gi 
(sec the r values in Tablc|T|. The neutron fit is completely 
dominated by the high-precision EOl-012 data, which are 
concentrated on a very limited Xb range; it clearly indi- 
cates a 18-40% breaking of the WW relation at high xb, 
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proton 


x7d-o.f. 


not 


^*low 




(I) A = 


1.22 








(II) A = a(l - xb) 13 ((/3 + 2)x B - 1) 
a = 0.13 ±0.05 










/3 = 4.4 ± 1.0 


1.05 


15-32% 


18-36% 


14-31% 


neutron 


(I) A = 


1.66 








(II) A = q(1 - xb) 13 (W ± 2)zb - 1) 
a = 0.64 ±0.92 










(3 = 24 ± 10 


1.11 






18-40% 



TABLE I: Results of the 1-parameter fits of the WW breaking term A for different choices of its functional form. The value r 
of the relative size of the breaking term is computed for three regions of xb'- the entire measured xb range, [0.02,1]; the \ow-xb 
region, [0.02,0.15]; and the high-XB region, [0.15,1]. See text for further details. 




FIG. 1: Top panels: Experimental proton and neutron <?2 structure functions compared to g? ■ The crosses represent <?2 
computed at the experimental kinematics, while the solid lines are gY W computed at the average Q 2 of the E155x experiment. 
Data points for the proton target [s| [4{| have been slightly shifted in xb for clarity. For the neutron only the high-precision 
data from Isol IHll ] are included. Bottom panels: The WW-breaking term A fitted to A cx computed using the LSS2006 gV 1 
(hashed region). The dashed line represents (?2 ~ (ffa )', the spurious HT contribution to A that would be obtained using 
the total gi to compute A cx . 



but cannot be used to conclude much at lower xb values. 
A striking feature of the proton WW-breaking term in 
Fig. [1] is that it is comparable in size and opposite in sign 
to g^ w — (32 VW )'- It is essential, therefore, to use fits 
of gi that subtract higher twist terms, which would oth- 
erwise largely cancel the proton WW-breaking term and 
obscure the violation of the WW relation. In the case of 
the neutron one would generally obtain an enhancement 
of the WW-breaking term, although the experimental un- 
certainties there are considerably larger. 

In summary, we have found that the experimental data 
are consistent with a substantial breaking of the WW 
relation |2]). Previous analyses have verified the WW 
relation only qualitatively, and using parametrizations 



which do not subtract higher twist terms in g\. The 
present analysis clearly demonstrates that this can give 
the misleading impression that the WW relation holds 
to much better accuracy than it does in more complete 
analyses where the higher twist corrections have been 
consistently taken into account. More data are certainly 
needed to pin down the breaking of the WW relation to 
higher precision. New data are expected soon from the 
HERMES Collaboration and from the d2n (E06-014)and 
SANE (E07-003) experiments at Jefferson Lab [H, ^1]. 
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IV. TOWARD A DEEPER UNDERSTANDING 
OF QUARK-GLUON-QUARK CORRELATIONS 

In the past, since the LIR-breaking gx term was not 
considered in Eq. (I4T[) and the quark-mass term with 
hi was neglected, the breaking of the WW relation was 
considered to be a direct measurement of the pure twist- 
3 term gx- The presumed experimental validity of the 
WW relation was therefore taken as evidence that gx is 
small. This observation was then generalized to assume 
that all pure twist-3 terms are small. In contrast, the 
present analysis shows that, precisely due to the pres- 
ence of gx, the measurement of the breaking of the WW 
relation does not provide information on a single pure 
twist-3 matrix element. Even if in future the WW rela- 
tion were to be found to be satisfied to greater accuracy 
than the present data suggest, one could only conclude 
that the sum of the terms in (j47|) is small, 

E e a(-3T^)+ / 1 -(?t(2/) + ^(2/))) «0. (54) 
a V Jx V J 

This can occur either because and g T are both small, 
or because they (accidentally) cancel each other. No 
information can be obtained on the size of the twist-3 
quark-gluon-quark term gx from the experimental data 
on <?2 alone. Note that these results were essentially al- 
ready obtained in Ref. [34]. In that work, however, the 
authors considered the WW breaking to be small and 
assumed that g T was small (which we argue is not neces- 
sarily the case), concluding that g T is also small. 

Of course it is desirable to test our conclusions em- 
pirically. A reliable way to investigate gx experimen- 
tally is through measurement of the function g[x- This 
function is accessible in semi-inclusive deep inelastic scat- 
tering with transversely polarized targets and longitudi- 
nally polarized lepton beams (see, e.g., the second line 
of Tab. IV in Ref. [Ilj]). Preliminary data related to this 
function have been presented by the COMPASS Collab- 
oration [5?! and more are expected from the HERMES 
Collaboration and from the E06-010 experiment at Jeffer- 
son Lab [58| . Using the EOM relation (|36| and assuming 
to = 0, one obtains 

xg T (x) = xg T {x) - g$(x) . (55) 

In combination with the measurement of the WW break- 
ing, this can be used to determine the size of twist-3 
function cjx- (Alternatively, one can use the LIR (|25p.) 

The importance of separately studying gx and <?t re- 
sides in the fact that these are projections of differ- 
ent combinations of the twist-3 functions Gf(x,x') and 
Gf(x,x'). As with all other terms in the decomposition 
of the quark-gluon-quark correlator in Eq. ([25]) . these 
functions are involved in the evolution equation of twist- 
3 collinear PDFs [5^, [6(3] , in the evolution of the trans- 
verse moments of the TMDs [6l|, [62j , in the calculation 
of processes at high transverse momentum [38| . and in 



the calculation of the high transverse momentum tails of 
TMDs [H [H. Ultimately, through a global study of 
all of these observables, one could simultaneously obtain 
better knowledge of twist-3 collinear functions and twist- 
2 TMDs, and at the same time test the validity of the 
formalism. Gathering as much information as one can 
on the quark-gluon-quark correlator is essential to reach 
this goal. The separation of the functions gx and 'gx is 
an important first step in this direction. 



V. CONCLUSIONS 

In this analysis we have shown that the Wandzura- 
Wilczek relation for the gi structure function is violated 
by a quark mass term, and two distinct pure twist-3 con- 
tributions, containing the parton distribution functions 
gx and gx- As evident from their definitions in Eqs. (|3ip 
and ([57]) respectively, these correspond to two different 
projections of the general quark-gluon-quark correlator 
in Eq. (|27p . Their measurement can give unique and 
complementary information on twist-3 physics. 

The two twist-3 functions have some interesting con- 
nections with the formalism of transverse momentum dis- 
tributions. One of them is involved in the equation-of- 
motion relation expressed in Eq. (|3l)|) . while the other is 
involved in the Lorentz invariance relation in Eq. (|25[) . 
Both relations contain the same moment of the trans- 
verse momentum distribution g\x- From the theoretical 
point of view, this is another intriguing example of the in- 
terplay between transverse momentum distributions and 
(collinear) twist-3 distributions. From the phenomcno- 
logical point of view, this means that a measurement of 
the function gyx in semi-inclusive DIS in principle allows 
one to separately measure gx and gx- 

Although the Wandzura-Wilczek relation is often used 
to simplify the treatment of twist-3 and TMD physics, 
we stress that there are no compelling theoretical or phe- 
nomenological grounds supporting its validity. In fact, 
using the experimental information currently available, 
we were able to provide a quantitative assessment of the 
violation of the Wandzura-Wilczek relation. Assuming 
a simple functional form for the WW-breaking term, we 
found that it can be as large as 15-40% at the 1-a con- 
fidence level. 

As new data become available, it should be possible to 
better pin down the violation of the Wandzura-Wilczek 
relation and measure the transverse momentum distri- 
bution gix in semi-inclusive DIS. This will offer us a 
deeper look into the physics of quark-gluon-quark corre- 
lations and its connection to transverse momentum dis- 
tributions. 
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APPENDIX A: TMDS WITH A NON-LIGHTLIKE 
WILSON LINE DIRECTION 



Factorization theorems beyond tree-level [2J, |26j, [27J, 
IH,[6^] demand a slightly non-lightlike vector v in order to 
regularize the lightcone (or rapidity) divergences [H[ . 
In Ref. [24] the Wilson line vector is chosen to be timelike 
and a parameter ( 2 = 4(P • v) 2 /v 2 is used as a regulator, 
with the requirement that £ 2 3> M 2 ,k\. In other articles 
in the literature v has been chosen to be spacelike 25] . 

In addition to k ■ P, k 2 , P ■ v and k ■ v, the PCFs A; 
and Bi can now in principle depend also on v 2 . We can 
derive the following relation between the invariants 



becomes 



k ■ v 
P ■ v 



2(7 



( 2 (l + a) ' 



(Al) 



with a = \J\ — 4M 2 /( 2 . Neglecting terms of order 
M 2 /£ 2 and cr/£ 2 , the above expression reduces to x. We 
therefore conclude that the PCFs depend on cr, r, x and 
additionally on £ 2 . To be precise, the definition of par- 
ton correlation functions in [27} involves an additional 
soft factor which is not included in the correlator $. The 
inclusion of the soft factor leads to an additional depen- 
dence on a gluon rapidity parameter. However, we leave 
this soft factor aside since it plays no role in our subse- 
quent discussion. 

The expressions for the TMDs in Eqs. (HU), ^ and 
then become 



g 1L (x,k 2 ,,C 2 ) = J dadr S(t - xa + x 2 M 2 + k 2 ,) 



A 7 + xA s 



AM 2 



C 2 (i + «) 



-Pi 



a - 2xM 2 
2M 2 



4M 2 



C 2 (i + «) 



B 



13 



(A2) 



gi T (x, k\, C 2 ) = J dadr S(t - xa + x 2 M 2 + fc 2 .) 



A 7 + xA s 



AM 2 



-B 



(A3) 



C 2 (l + a) 

g T (x, fef,, C 2 ) = J dadr 5(t - xa + x 2 M 2 + k\) 

. (A4) 



t — xa + x 2 M 2 . 
~ A 6 ^779 ^8 



2M' 



The full expression for c/t which generalizes Eq. (|26[) then 



g T (x) = J d 2 k T dadr S(t - xa + x 2 M 2 + fcf,) 



4M 2 

B n + xB 12 + — B 14 
C 2 (l + a) 



fc| /dA 7 
~ 2M 2 V dx 



dA 8 

■ x_ 5 1" 

ox 



AM 2 dB 



13 



C 2 (l + a) dx 



7T J dadr S(t - xa + x 2 M 2 + fe|) fe| 



a - 2xM 2 ( . 



4M 2 



C 2 (l + a) 



Bis) 



k T — >oo 

(A5) 



APPENDIX B: PARTON CORRELATION 
FUNCTIONS FOR A QUARK TARGET 

In this Appendix we compute the parton correlation 
functions relevant for our discussion of the WW relation 
for the case of a point-like quark target. The calculations 
are performed in the first non-trivial order in perturba- 
tive QCD (i.e., at order a s ) [13, H]. To this end we 
insert a complete set of intermediate states into Eq. ((U) . 
To order a s , only the vacuum state and a one-gluon state 
are relevant. The involved Feynman diagrams are shown 
in Fig. [2] (real gluon contributions) and Fig. [3] (virtual 
gluon contributions). 

The correlator may be written as 



P, S; v) = S^(P - fc)0>™(™ 2 , A 2 , C 2 , A) 



$lf(k,P,S;v), 



(Bl) 



where $ vlr denotes the contributions from the vacuum 
intermediate state. Its kinematics is totally determined 
by the four-dimensional delta-function 5^ (P — k) and 
depends only on the quark mass m, with a small gluon 
mass A serving here as an infrared regulator, and the 
parameter ( 2 = 4(P • v) 2 /v 2 which regulates lightcone 
divergences. By applying a renormalization procedure 
we can subtract ultra-violet divergences in <!> vir , which 
introduces a dependence on the renormalization point 
H R . The virtual corrections can be written as 

$T/(fe, P, S; v) = S^(P - k)(P, S, d\ ^(0) w ( " 0iOO) |0) 
x (OlW^o) ^(0)\P,S,d), 

(B2) 

where the incoming on-shell quark is described by the 
state |P, S, d) , with d a color index of the quark in the fun- 
damental SU(3) representation. For the sake of brevity 
we will omit the explicit dependence on and summation 
over the color indices in the following. Since we work in 
Feynman gauge, possible contributions from gauge links 
at lightcone infinity are irrelevant [L6| . 

The second contribution in Eq. (|B1[) is generated by 
one gluon in the intermediate state. To order a s it is 
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TTG1S> 




FIG. 2: Diagrams in the quark-target calculation involving only real gluons. The Hermitean conjugate diagrams, which are 
not shown, are also taken into account in the calculation. 




FIG. 3: As in Fig. [2] but for diagrams involving virtual gluons. 



given by 



$lf(k, P, S; v) = 5+ (( p k ) 2 A ') 



X M"/ H (k, P, S; v) M°' P {k, P, S; v) 



(B3) 



with M = Mt 7 °, S+(a 2 ) ee 6(a 2 )6(a°), a denotes the 
polarization of the gluon in the intermediate state, and 
(3 is its color index in the adjoint representation of SU(3). 
The matrix element M is then represented by 



M?' (k, P, S; v) = (P- k, a, /3|^(0)|P, S, d) 

+ ig dX{P-k,a,P\v-A(Xv)i/;i(0)\P,S,d), 
Jo 

(B4) 



where \P — k, cr, f3) denotes the intermediate gluon state 
with a color index (3. The leading perturbative contribu- 
tion in a s to the matrix element M gives 



M?' f, (k,P,S;v) = -gif 



[k 2 



m z + ie\ 



ve* a (P-k) 
[v(P-k)+ie]J u 



ui(P,S), 



(B5) 



where e(P — k) denotes the gluon polarization vector and 
u is the quark spinor. The color flow is given by the 
color matrix t@ in the fundamental representation. In- 



serting (|B5|) into (|B3|) then yields 
$lf(k, P, S; v) = -^C F S+((P - k) 2 A 2 ) 

{ft + m) 7M (f + m)(l + 75^)7^ + m) 



[k 2 



e][k 2 



(f + m)(l + 75 ^(^ + m) 



[k 2 



ie][v -(P-k) + if] 



(ft+m^(f + m)(l + 75 ft) 
[k 2 — m 2 + ie] [v ■ (P — k) — ie] 

v 2 (f + m)(l+ l5 $) 



[v ■ (P - k) +ie][v (P - k) 



ie 



(B6) 



The various parton correlation functions in Eq. ([5]) can be 
extracted from Eq. (;BGj by decomposing the numerators 
in terms of the basis matrices 1, 75, 7 M , 7^75 and a^ u . 
In this way we obtain expressions for parton correlation 
functions at leading order in a s for a quark target. In the 
following we list only the PCFs ^6-8 and Pn_i4 which 
are relevant for the discussion of the Wandz ura-Wilczek 
relation, cf. Eqs. flH])-([II])- Setting a = Vl - 4to 2 /C 2 , 
we find (to order a s ) 

CFOis 



A%*(t, a, x, C 2 ) = ^S+(r -a + m 2 -X 2 ) 

T + m 2 (1 + a)(l + ax) + 2a/( 2 

r - m 2 ) 2 + [r- ™ 2 ] [(1 + a)(l - ax) - 2a /( 2 ] 

2(1 + a) 2 

[(1 - a.x) 2 (l + a) 2 C 2 - ia(l - ax)(l + a) + 4a 2 /( 2 ] 

(B7) 



A*r\r,a,x,( 2 ) = 0, 
~2m 2 



5+(i 



A 2 ) 



(B8) 



(B9) 



11 



Cfchi 

~2^ 



5+(t - a + m 2 - A 2 ) 
-(1 + a) 



[r - m 2 ] [(1 + o)(l - ax) - 2ct/C 2 



2tt 2 



o+(t-<t + to 2 - A 2 ) 



(1 + a) 



[r - m 2 ] [(1 + o)(l - ax) - 2<j/C 2 } 
B[T(T,a,xX 



' - ~ " • ' ~~ [T — a + m — A ) 



2tt 2 



-(1 + a) 



(BIO) 
(Bll) 

(B12) 
(B13) 



[r - to 2 ] [(1 + a)(l - ax) - 2a/C 2 ] 
^f 1 (r,(7,a ; ,C 2 )=0. 



These results demonstrate that all terms in Eq. (|A5[) con- 
tribute to generate a nonzero gr since (i) the -B^ terms are 
nonzero, (ii) the PCFs can depend explicitly on x, and 
(iii) the boundary term at k\ = oo cannot be neglected. 



APPENDIX C: QUARK TARGET TMDS AND 
PDFS AT x < 1 

We are now in a position to calculate the TMDs for 
a quark target defined in Eqs. (|A2|) (|A4|) . their fc^- 
integrals appearing in the LIR of Eq. (123)) , and the func- 
tion gx as defined in Eq. (|A5| . Similar calculations have 
been performed in [24], l37ll44l l69l [70j . Without entering 
into details, we note that the light-cone divergences oc- 
curring for £ — > oo can be moved to x = 1, introducing 
the well-known "plus" distribution [H, [33J . If we restrict 
ourselves to the region x < 1, the results are free of light- 
cone divergences and do not depend on £. In this region 
we can use either Eqs. (TA"2"|) - (rM|) or — (PJ3J) - The re- 
sulting functions are then given by 



g 1L (x < 1,*4) 



2C F a s 



1 



(27r) 2 k\ + xX 2 + (1 - x) 2 m 2 
2(1 - x)(l - x(l - x))m 2 2x 



fc 2 . + rrA 2 + (1 - x) 2 m 2 



(1 



(CI) 



ai T (a; < l,fefO 



2Cfck s 



2x(l — x)m 2 



g T (x<\,k 2 T) 



(2tt) 2 (fc 2 + xX 2 + (1 - x) 2 m 2 ) 2 ' 

(C2) 

1 



(2tt) 2 k 2 + x\ 2 + (1 - .X) 2 77. 
(1 - x) 2 (l + x)m 2 1 + x 

k\ + xX 2 + (1 - x) 2 m 2 (l-x) + 



(C3) 



When working with non-lightlike Wilson lines, it is not 
clear how to obtain the collinear parton distribution 
functions upon integration over the transverse momen- 
tum [24j. However, at the one-loop level these subtleties 
are relevant only at x = 1. Since we restrict ourselves to 
the region x < 1, we can safely compute collinear PDFs 
through fc^-integration. For simplicity we choose an up- 
per boundary Q for the fey-integration, and shift quark 
mass effects into the finite part by introducing an arbi- 
trary infrared cutoff i£_jn order to obtain agreement with 
the results of Refs. [37], The divergent parts of the 
parton distributions, i.e., the terms including the upper 
cutoff Q, are given by 



. a s C F 1 + x 2 Q 2 
g 1L (x < 1 = — — — In — , 

2tt (1 — x) + [i £ 

. a s C F l + 2x-x 2 O 2 
g T (x < 1) = — ; In- 



2tt (1~x)- 



9it( x < 1 ) = ~~lT~r x ( 1 ~ x)ln^j. 

A* 



2tt 



(C4) 
(C5) 
(C6) 



These results have appeared earlier in Refs. [2J, [33, 
HI, Iz3|, but have been derived here for the first time 
starting from the PCFs. 

For ~gi at x < 1, using either Eq. (|A5|) or Eq. (f2l)|) we 
obtain 



~ , . ^ asC F Q 2 
g T (x < 1 = — — In — , 
2n n £ 



(C7) 



confirming the result in Eq. (|35[) . which was not obtained 
directly but rather using the LIR relation Eq. ([2"5]) . 
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